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Free energy of a static quark anti-quark pair and the renormalized Polyakov loop in 
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We study the free energy of a static quark anti-quark (QQ) pair at finite temperature in three 
flavor QCD with degenerate quark masses using N T = 4 and 6 lattices with Asqtad staggered fermion 
action. The static free energy was calculated for different values of the quark mass and the entropy 
contribution at large distances has been extracted. We also calculate the renormalized Polyakov 
loop following the approach by Kaczmarek et al. 



I. INTRODUCTION 

o ■ ...... . r 

Non-perturbativcly the in-medium modification of inter-quark forces, e.g. screenirigof fundamental charges at finite 
£NJ | temperature, is usually studied in terms of the free energy of static (anti)quarks 0, UIE QL IE IE • More precisely one 
calculates the difference in the free energy of the system with static quarks and the same system without quarks at 
fixed temperature. The free energy of static quark anti-quark for different color channels (singlet, octet and averaged) 
has been studied in great detail in quenched QCD (pure SU(3) gauge theory) 0,0,0 as well as m SU(2) gauge theory 
(E0- In the case of full QCD only the color averaged free energy was extensively studied; the first results for singlet 
and octet free energy for two flavor QCD have appeared only very recently 0. In this paper we are going to study 
the free energy of static quark anti-quark pair in 3 flavor QCD using the so-called Asqtad staggered fermion action 
! with two different lattice spacings (corresponding to Nt — 4 and 6) at three different quark masses. 

Apart from the in-medium modification of inter-quark forces the study of the static free energy is interesting as it 
Q\ \ can be used to define the renormalized Polyakov loop. Although the Polyakov loop is not an order parameter in the 
^) • presence of dynamical quarks it shows rapid variation in the transition region and therefore is widely used to describe 
the transition (crossover) in full QCD. In particular, they are useful for constructing effective mean-field theories 
and studying the interplay between the chiral and deconfining aspects of the transition 10|. We will study the 
temperature dependence of the renormalized Polyakov loop which shares most of the properties of the usual Polyakov 
loop but has a meaningful continuum limit. 

The rest of the paper is organized as follows. In section 11 we discuss the lattice setup, parameters of simulations 
and the zero temperature potential which turns out to be crucial for the analysis of the finite temperature free energies 
'"" j" 1 ■ presented in section III. Conclusions are presented in section IV. 

^ ■ II. PARAMETERS OF SIMULATIONS AND THE STATIC POTENTIAL AT ZERO TEMPERATURE 
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In this section we are going to describe the lattice parameters for which our analysis have been performed. In 
our study we use staggered fermions with Asqtad action 0. Our analysis to a large extent is based on the gauge 
configurations generated by the MILC collaboration usin g th e Asqtad action. Therefore we adopt their strategy for 
fixing the parameters which is described in Refs. [Ill Il3j . namely the strange quark mass was fixed by requiring 
that m^/m^ = 0.673, and the temperature scale (the inverse lattice spacing) was fixed from the scale n defined from 
the zero temperature static potential as 

r2 ^L, - L (i) 

The RG inspired Ansatz for the gauge coupling dependence of the lattice spacing Eil was used for the scale setting. 
We use the most recent value of r\ extrapolated to continuum and to the physical value of the light quark masses 
r*i = 0.317 fm 15] to convert the lattice units to physical units. In Table[I]we summarize the lattices, quark masses 
and the corresponding temperature range used in our analysis. 

The free energy of a static quark anti-quark pair contains a lattice spacing dependent divergent piece and thus 
needs to be renormalized. This can be done by normalizing it to the zero temperature potential at short distances 
where the temperature dependence of the free energy can be neglected 0. The static quark potential has been studied 
by the MILC collaboration at three different lattice spacings and various quark masses 
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The zero 

temperature potential is also defined only up to some constant which needs to be fixed. As temperature is varied by 
varying the lattice spacing we need to specify the form of the zero temperature potential which is valid for the whole 
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FIG. 1: The static potential at zero temperature from the MILC collaboration versus the parameterization @ (right) and the 
effective coupling constant a a (r) as a function of distance (left). Both cases of three degenerate quarks denoted as 3F as well 
as one heavy (m q = m 3 ) and two light denoted as 2+1F are shown. 



lattice 


miight 


T 


12 3 x 4 


0.2m s , 0Am s , 0.6m s 


135 - 393 MeV 


8 3 x 4 


0.2m s , 0.4m s , 0.6m s 


135 - 412 MeV 


12 3 x 6 


0.2m s , 0.4m s , 0.6m s 


145 - 310MeV 



TABLE I: The parameters of simulations, m s denotes the strange quark mass. 



range of lattice spacing relevant in this study, i.e. 0.09/m < a < 0.3/to. We choose the following form of the zero 
temperature potential 

O- 44 „ „ 0.0125 
nV(x) = h 0.56 ■ x H rj— , x = r/n (2) 

X X 

In Fig. we show this parametrization of the potential against lattice data at three different lattice spacing, and 
different quark masses. Although in the present paper we discuss the case of three degenerate flavors, we show the 
zero temperature potential also for 2 + 1 flavor case as here more detailed data are available. Because of the very weak 
quark mass dependence the case of 2 + 1 flavors can be used as a good reference. As one can see this parameterization 
gives a fare description of the data for all lattice spacings and quark masses 01 • The last term in Eq. mimics the 
effect of the running coupling. In Fig. ^we also show the effective coupling constant a s (r) defined as 

3 2 dV(r) 

One can see again that Ansatz Q gives a good description of the lattice data. We note that in order to reduce the 
effects of the breaking of the rotational symmetry in the heavy quark potential following Refs. 0, the separation 
r between the static charges was redefined as r = n — (47rC , L(r))~ 1 , where Ci(r) is the lattice Coulomb potential 
for the Luescher-Weisz action In what follows r will always refer to the modified separation defined above. 



III. THE FREE ENERGY OF STATIC QUARK ANTI-QUARK 



Following Refs. |2Cj, J2JJ| the free energy of static quark- antiquark ( QQ) pair, i.e. the free energy difference of the 
system with and without static sources at fixed temperature T in the color singlet and octet channels is defined as 

eM-Fi (r, T)/T + C) = ^Tv(W(r)W^ (0)) (4) 

exp(-F 8 (r,T)/T + C) = -(TW^TW^O)) - — Tr (W (r)W^ (G)) (5) 

8 24 

Here W(x) = nfio 1 Uq(t,x) is the temporal Wilson line and L(x) — TrW(x) is known as the Polyakov loop. As 
W{x) is not gauge invariant one needs to fix a gauge in order to calculate F\ and Fg. In this study we use the 
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FIG. 2: The singlet (left) and octet (right) free energy of static QQ pair calculated on 12 3 x 4 lattice at m q = 0.4m s at different 
temperature in MeV. The solid line is the parameterization of the zero temperature potential by Eq. [5] 



Coulomb gauge as for this gauge a transfer matrix can be constructed and in the zero temperature limit the usual 
static potential will be recovered. Alternatively one can replace the Wilson line W(x) by a dressed gauge invariant 
Wilson line W(x) using the eigenvectors of covariant spatial Laplacian [3j. The dressed Wilson line, however, is a 
non-local operator. Both definitions turned out to equivalent Q (at least numerically). One can also consider the 
color averaged free energy defined as 

cM-Fav(r,T)/T + C) = i(TW(r)TWt( )) = ^{L(r)L\0)} (6) 

which can be written as a thermal average of the singlet and octet free energies 

eM-F av (r,T)/T) = ^exp{-F 1 (r,T)/T) + | cxp(-F 8 (r, T)/T). (7) 

This quantity is expressed in terms of local explicitly gauge invariant operators and for this reason was extensively 
studied in the past also in full QCD [2^, [23| . The normalization constant C needs to be fixed using some physical 
normalization condition. As we expect that at very short distances the free energy of static quark anti-quark pair is 
only given by their interaction energy in the vacuum we fix C by requiring that the singlet free energy coincides with 
the zero temperature potential given by Eq. J2J). 

We start the discussion of our numerical results with the case of the quark of mass 0.4m s on the 12 3 x 4 lattice [27) . 
The corresponding numerical results for the singlet and octet free energy are shown in Fig. [21 The singlet free energy 
approaches a finite value at large distances which is usually interpreted as string breaking at low temperature and 
screening at high ones. Note that the distance where the free energy effectively flattens is temperature dependent; it 
becomes smaller at higher temperatures. At small distances the singlet free energy is temperature independent and 
coincides with the zero temperature potential. This is intuitively expected as at small distances medium effects are 
not important. Similar observation has been made for quenched QCD |4J, |y,l7j. The octet free energy shows much 
stronger temperature dependence. At short distances it is expected to have a repulsive tail, which is clearly visible at 
high temperatures. At low temperatures the presence of such a repulsive tail is less obvious. The reason for this is 
the following. At low temperatures no data is available at really short distances (< 0.2/m) and the (large distance ) 
asymptotic value of the octet free energy is large compared to the value of the repulsive tail (which naively is a s / (6r)). 
Also the statistical accuracy is lowest at low temperatures. As the temperature increases both the value of the free 
energy at large distances is smaller and the statistical accuracy is higher; in addition more data at shorter distances 
become accessible. The numerical results for other values of the quark masses are similar. 

As N T — 4 lattices correspond to a quite large lattice spacing, one may worry about lattice artifacts in the static 
free energies. Therefore we have compared the singlet free energies calculated from N T = 4 and N T — 6 lattices at 
approximately same temperatures as shown in Fig. [3J As one can see, there is no sizeable lattice spacing dependence 
even at the lowest temperatures where the lattice spacing is the largest. Thus calculations of the quark anti-quark 
free energy on N T = 4 lattices are justified. 

We also investigate the temperature dependence of the color averaged free energy. In general due to the presence 
of color octet contribution the temperature dependence of the color averaged free energy is stronger than that of the 
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FIG. 3: The singlet free energy calculated on 8 3 x 4, 12 3 x 4 and 12 3 x 6 lattices at three different temperatures. 
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FIG. 4: The color averaged free energy at different temperatures calculated on N T = 4 (open symbols) and N T = 6 (filled 
symbols) lattices. The solid line is the zero temperature potential. 



singlet free energy |23|. At small distances, however, it is expected to be dominated by color singlet contribution and 
the approximate relation F av (r, T) ~ F t (r, T) + T In 9 should hold 0] . Therefore in Fig. H we show F av (r,T)-T In 9 
at different temperatures including results both from N T = 4 and 6 lattices. As one can see from the figure for 
temperatures T < TJAMeV F av — Tin 9 coincides with the zero temperature potential at the shortest distance 
available in this study. As at the shortest distance the zero temperature potential coincides with the singlet free 
energy (due to normalization) this in turn implies that the relation F av — Fx + T In 9 holds at this distance. 

As expected, the temperature dependence of the color averaged free energy is stronger than of the color singlet one. 

From Figs l2l4l we see that the free energy for all color channels reach a constant value at large distances F^T) = 
lim r _ 00 Fi(r, T), % = 1, 8, av. In the infinite volume limit we expect that F^(T) = F^T) — F™(T) as at very large 
distances the free energy of a static quark and anti-quark should not depend on their relative color orientation. This 
is in fact confirmed in the deconfined phase of pure gauge theory On our small lattices, 8 3 x 4, 12 3 x 6 we 

see small differences between F^ and F^ av which vanish within statistical errors when we go to 12 3 x 4 lattice. In 
practice we determine F^iT) by fitting the corresponding free energies by a constant at large distances. The lower 
limit of the fit range was always chosen such that the determined value of F^ does not depend on it within our 
statistical accuracy. 

The Polyakov loop L(x) — TrW(x) is the order parameter for deconfining transition in pure gauge theories. In 
full QCD, where dynamical quarks are present, it is no longer an order parameter as it has non-zero value in low 
temperature phase. Nevertheless, it is used to study the deconfining aspects of the transition in full QCD as well in 
effective models 0, ^| . The correlator of Polyakov loops corresponding to the color averaged free energy satisfies the 
cluster decomposition 

(L(r)Lt( ))|_ oo = |(L(0))| 2 . (8) 
Therefore following Ref . , we can define the renormalized Polyakov loop as 

F av (T) 

L ren = exp(-^i). (9) 
The numerical results for L ren (T) will be discussed in detail at the end of this section. In three flavor QCD for 
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FIG. 5: The effective screening radius from the singlet free energy as a function of the temperature for three different quark 
masses (left) and as function of T/T c for 0.4m s and quenched QCD (right) 



the quark masses used in this study there is no phase transition but only a crossover |ll|, [12|, [13|, |24| . Nevertheless, 
both the chiral condensate and the unrenormalized Polyakov loop show a rapid change at approximately the same 
temperature [Til IrH Il3l l23j referred to as the transition temperature T c . Therefore we can define the transition 
temperature as the temperature where dL ren (T)/dT has a maximum, or equivalently dF ren (T) / dT has a minimum. 
We have determined T c in this way and have found it to be consistent with the value of T c defined from the peak of 
the chiral susceptibility jllL Hia [jjj . In the following discussion we will find instructive to plot different quantities as 
function of T/T c , especially when comparing with calculations done in SU(3) gauge theory. 

To characterize the range of interaction in the medium it is convenient to introduce the effective screening radius 
r scr . It is defined as the distance at which the singlet free energy is only 10% below its asymptotic value F\{r = 
r S cr,T) — 0.9F^ o (T). Here F^(T) is the asymptotic value of the singlet free energy at infinite separation. In Fig. [3] 
we show the values of r scr for three different quark masses and 12 3 x 4 lattices. Certainly as F\(r,T) has statistical 
errors, it is difficult to determine exactly at which distance r the equation F\(r = r scr ,T) — Q.9F^(T) holds. We have 
tried to estimate this uncertainties in the values of r scr and show them in Fig. [SJas errors bars. At small temperatures 
the value of the screening radius is about 0.9 fm and is temperature independent. As we increase the temperature 
r scr decreases reaching the value of 0.5 fm at the highest temperature. Note that the temperature dependence of r scr 
is roughly the same for all quark masses. In Fig. [SJwe also show the effective screening radius for 0.4m s and quenched 
QCD estimated from the singlet free energies of Ref. 4] as function of T/T c . Close to T c the screening radius in 
quenched QCD is considerably higher than in three flavor QCD but at higher temperatures they are comparable. 

Now we are going to investigate the temperature and quark mass dependence of the asymptotic value of the free 
energy and identify the entropy contribution to it. In the following discussion we will take the asymptotic value of 
the free energy in the infinite volume limit to be given by the color averaged value F^ (T) = F^ (T) and also skip the 
index i in the following discussion. The numerical results for i 7 ^ (T) at three different quark masses are shown in Fig. 
[S] One immediately notices that the quark mass dependence may vanish at small temperatures (T < 150AleV) and 
definitely negligible at high temperatures (T > 250MeV) while large mass dependence is observed in the transition 
region. One may wonder to which extent this mass dependence is due to the mass dependence of the transition 
temperature T c . Therefore in Fig. |H|we also show the data for F^T) versus T/T c which shows that up to T/T c — 1 
there is no mass dependence. This is in accordance with finding of Ref. |23| where for quark masses below the strange 
quark mass the mass dependence of F^ is quite small for T/T c close to one. Thus for temperatures below T c the 
mass dependence of Foo(T) can be understood in terms of mass dependence of T c . However, as one can also see from 
Fig. El for T/T c > 1 substantial mass dependence is seen. 

In the case of very small temperatures we expect Foo(T) to be temperature independent and related to twice the 
binding energy of a heavy-light (D— or B-) meson 2i?^„ = 2Md,b — 2m c ^. More precisely, it should be the binding 
energy of a static-light system as found in the 3d SU(2) Higgs model which is quite similar to QCD [25|. Because of 
the heavy quark symmetry we expect the binding energy of a heavy-light meson to be the same as of the static-light 
system. Based on this observation it has been argued in Ref. 26] that the decrease of Foo(T) with the temperature 
close to T c implies the decrease of the M^^b leading to quarkonium suppression. However, one should keep in mind 
that Foo(T) also contains an entropy contribution. The entropy due to the presence of a static QQ pair can be 
calculated as 
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FIG. 6: The asymptotic value of the free energy at large distance for different quark masses as function of the temrjerature in 
physical units(left) and of T/T c (right). On the right plot we also show the value of Foo from pure gauge theory 
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FIG. 7: The asymptotic values of the entropy and the internal energy of static QQ pair at large distances as function of the 
temperature. 



We also can calculate the energy induced by a static QQ pair Uoo — Foo + TSoo. Numerically the derivative with 
respect to the temperature in Eq. i fTUj l was estimated using forward differences. In Fig. Qwe show the entropy 
and energy Uoo as functions of temperature. Both the entropy and the energy show a strong increase near T c and 
with our definition of T c peak exactly at T c . This large increase in entropy and energy is probably due to many-body 
effects and makes the interpretation of Uoo as the binding energy of heavy-light meson not very plausible. 

The numerical results for Foo(T) can be easily converted to the numerical results for the renormalized Polyakov 
loop L ren (T) which are shown in Fig. [SJ The mass dependence of L ren (T) is the same as for Foe, when plotted versus 
temperature in physical units it shows mass dependence for lbQMeV < T < 250MeV~ and when plotted as function 
of T/T c mass dependence for T/T c > 1. The renormalized Polyakov loop should not depend on the lattices spacing. 
To check this we have calculated L ren (T) also on 12 3 x 6 and 8 3 x 4 lattices. The results are summarized in Fig. [SJ 
and as one can see the N T (lattice spacing) dependence of L ren (T) is quite small. 



IV. CONCLUSIONS 



In this paper we have presented calculations of the free energy of a static quark anti-quark pair in three flavor QCD 
for several quark masses. We have found that the free energy gets screened beyond some distance for all temperatures, 
as expected. For small temperature this distance, the effective screening radius, does not depend on the temperature 
and is about 0.9 fm. As the temperature increases the effective screening radius decreases. Using the asymptotic value 
of the free energy we have defined the renormalized Polyakov loop and proved its scaling with the lattice spacing. 
We have also identified the entropy contribution to the free energy as well as the internal energy at large distances 
and found that they show strong increase at T c . We have found substantial quark mass dependence in the vicinity of 
the transition. In the future it will be certainly interesting to extend this study to larger lattices and to study the 
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FIG. 8: The renormalized Polyakov loop calculated on 12 3 x 4 lattices for three different quark masses (left) and for m q = 
0.2m s ,0.6m s on 8 3 x 4, 12 3 x 4 and 12 3 x 6 lattices (right). The insertion in the left figure shows the renormalized Polyakov 
loop as function of T/T c (see text). . 



dependence of the entropy contribution as function of the separation which because of the limited statistics was not 
possible here. Many properties of the static quark anti-quark free energies presented in this three flavor study turns 
out to be similar to the preliminary findings in two flavor case by the Bielefeld group reported in Ref. Q. 
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